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The classical Kramer sampling theorem establishes general conditions 
that allow the reconstruction of functions by mean of orthogonal sam- 
pling formulae. One major task in sampling theory is to find concrete, 
non trivial realizations of this theorem. In this paper we provide a new 
approach to this subject on the basis of the M. G. Krein's theory of 
representation of simple regular symmetric operators having deficiency 
indices (1, 1). We show that the resulting sampling formulae have the 
form of Lagrange interpolation series. We also characterize the space 
of functions reconstructible by our sampling formulae. Our construc- 
tion allows a rigorous treatment of certain ideas proposed recently in 
quantum gravity. 
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1 Introduction 



It has been argued recently that sampling theory might be the bridge that would allow to 
reconcile the continuous nature of physical fields with the need of discretization of space- 
time, as required by a yet-to-be- formulated theory of quantum gravity [TBI [HI EH]. This 
idea, which is partially developed in [T3] for the one-dimensional case, introduces the use 
of simple regular symmetric operators with deficiency indices (1,1) to obtain orthogonal 
sampling formulae. It is remarkable that the class of operators under consideration in [13] 
had been already studied in detail by M. G. Krein [HJ [TSl [20] more than 60 years ago. 

The conjunction of the ideas of [H] and Krein's theory of symmetric operators with 
equal deficiency indices suggests the means of treating sampling theory in a new and 
straightforward way. By introducing this new approach, we put in a mathematically rig- 
orous framework the ideas in [13] and propose a general method for obtaining analytic 
sampling formulae associated with the self-adjoint extensions of a simple regular symmet- 
ric operator with deficiency indices (1, 1). 

A seminal result in sampling theory is the Whittaker-Shannon-Kotel'nikov (WSK) sam- 
pling theorem [T6], [221 El] • This theorem states that functions that belong to a Paley- Wiener 
space may be uniquely reconstructed from their values at certain discrete sets of points. 
A general approach to WSK type formulae was developed by Kramer [17] based on the 
following result. Given a finite interval I = [a,b], let K(y,x) G L 2 (I,dy) for all x G R. 
Assume that there exists a sequence {x n }nez for which {K(y, x n )} n< zz forms an orthogonal 
basis of L 2 (I, dy). Let / be any function of the form 



for some g G L 2 (I, dy), where the brackets denotes the inner product in L 2 (I, dy). Then / 
can be reconstructed from its samples {f(x n )} ne z by the formula 



The search for concrete realizations of the Kramer sampling theorem, which in some 
cases has proven to be a difficult task, has motivated a large amount of literature. See, for 
instance, [31 El El HHJ [25j [26] and references therein. On the basis of the approach proposed 
in the present work, we obtain analytic sampling formulae for functions belonging to linear 
sets of analytic functions determined by a given regular simple symmetric operator. 

A central idea in Krein's theory is the fact that any simple symmetric operator with 
deficiency indices (1,1), acting in a certain Hilbert space TC, defines a bijective mapping 
from 7i onto a space of scalar functions of one complex variable having certain analytic 
properties. In this space of functions one can introduce an inner product. In the particular 
case when the starting point is an entire operator (see Definition 2), the space of functions 



f(x) := (K(.,x),g(.)) L2 
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turns out to be a Hilbert space of entire functions known as a de Branges space [I]. 

Summing up, in this paper we present an original technique in sampling theory based on 
Krein's theory of regular symmetric operators. Concretely, we obtain an analytic sampling 
formula valid for functions belonging to linear sets of analytic functions associated with 
symmetric operators of the kind already mentioned. We then focus our attention to entire 
operators to provide a characterization of the corresponding spaces of entire functions. As 
a byproduct, we also provide rigorous proofs to some of the results formally obtained in 
[T3j . In a sense this work may be considered as introductory to the theoretical framework. 
Further rigorous results related to [IB] , as well as other applications to sampling theory 
and inverse spectral theory, will be discussed in subsequent papers. 

Our approach to sampling theory allows us to reconstruct functions that are out of the 
scope of certain techniques developed recently [Sj. It is worth remarking that, in terms 
of the method derived in this paper, the results of [8] correspond to considering only a 
particular self-adjoint extension of an entire operator. We also remark that [8] does not 
fit in the framework proposed in [13] for which, in contrast, our results indeed give some 
rigorous justification. 

This paper is organized as follows. In Section 2, we introduce some concepts of the 
theory of regular symmetric operators. In Section 3, we state and prove a sampling formula. 
In Section 4, we provide a characterization of the Hilbert spaces associated with the class of 
operators under consideration, paying particular attention to the case of entire operators. 
Finally, in Section 5, we discuss some examples. 

2 Preliminaries 

Most of the mathematical background required in this work is based on Krein's theory of 
representation of symmetric operators as accounted in the expository book [TT| . In this 
section we review some definitions and results from the theory of symmetric operators and 
introduce the notation. 

Let TL denote a separable Hilbert space whose inner product (•, •) will be assumed 
anti-linear in its first argument. 

Definition 1. A closed symmetric operator A with domain and range in 7i is called: 

(a) simple if it does not have non-trivial invariant subspaces on which A is self-adjoint, 

(b) regular if every point in C is a point of regular type for A, that is, if the operator 
(A — zl) has bounded inverse for every zgC. 

By [TT1 Theorem 1.2.1], an operator A is simple if and only if 

p| Ran (A -zl) = {0}. 

z:Im z^O 
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Moreover, as shown in (TTJ Propositions 1.3.3, 1.3.5 and 1.3.6], a simple closed symmetric 
operator with deficiency indices (1, 1) is regular if and only if some (hence every) self-adjoint 
extension of A within 7i has only discrete spectrum of multiplicity 1. Also, 

dim [Ran {A - zl)] 1 - = dim [Ker (A* - zl)} = 1 (2.1) 

for all zgC For this kind of operators it is also known that, for any given igK, there is 
exactly one self-adjoint extension within 7i such that x is in its spectrum [TT| Proposition 
1.3.5]. 

Remark 1. In what follows, whenever we consider self-adjoint extensions of a symmetric 
operator A, we will always mean the self-adjoint restrictions of A*, in other words, the 
self-adjoint extensions of A within 7i (cf. Naimark's theory on generalized self-adjoint 
extensions [21 Appendix 1]). 

We will denote by Sym^' 1 \hi) the class of regular, simple, closed symmetric operators, 
defined on H, with deficiency indices (1,1). 

Let Aft be some self-adjoint extension of A G Sym^' 1 ^?-^). The generalized Cayley 
transform is defined as 

(At-wl) {A t -ziy\ 

for every w G C and z G C\ Sp(Aj). This operator has several properties [TU pag. 9]. We 
only mention the following one: 

(A t - wl) (Aft - zl)' 1 : Ker (A* - wl) -> Ker (A* - zl) (2.2) 

one-to-one and onto. 

A complex conjugation on 7i is a bijective anti-linear operator C : 7i — > TC such that 
C 2 = I and {Crj, dp) = (ip, rj) for all T], <p G TC. A symmetric operator A is said to be 
real with respect to a complex conjugation C if CT>om.(A) C Dom(A) and CAip = ACip 
for every ip G Dom(A). Clearly, the condition CDom(^4) C Dom(v4), along with C 2 = I, 
implies that CDom(A) = Dom(y4). If moreover A has deficiency indices (1,1), then A* is 
also real with respect to C (see the proof of (23J, Corollary 2.5]). 

3 Sampling formulae 

Let us consider an operator A G Sym^' 1 ' ) (7i). Let Aft be some self-adjoint extension of A. 
Given z G C \ Sp(Ajj) and ip G Ker (A* — z I), define 

V>(«) := {Aft - zqI) {Aft - zl)' 1 ^ = ^ + (z- z ) {Aft - zl)' 1 (3.1) 

for every z G C \ Sp(Ajj). This vector- valued function is analytic in the resolvent set of 
Aft and, by (12.21) . takes values in Ker {A* — zl) when evaluated at z. Moreover, ip{z) has 
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simple poles at the points of Sp(Afl) = Sp disc (/Ljj). By looking at the first equality in (13. ip . 
it is clear that the dependency of ip(z) on zq and i/jq is rather unessential. Indeed, for any 
other z' Q G C \ Sp(^4(j), we can take ip' = (Aft — z I) (Aft — z' I)~ ip in Ker (A* — z' I) and 
thus ip(z) = ip' Q + (z — z' Q ) (A$ — zl)" 1 ipQ by the first resolvent identity. 

Given z\ G C \ Sp(Ajj), let us choose some fi G TL such that (i[)(~zT), fi) ^ 0. The 
inner product (ip(z), jj) then defines an analytic function in C \ Sp(y4j), having zeroes at a 
countable set devoid of accumulation points in C. The set is defined as the subset 
of C for which 7i can not be written as the direct sum of Ran(A — zl) and Span{/i}. 
Following Krein, we call the element fi a gauge |2U] . 

Let us define _ 

fW - 0.2) 

(H,ip(z)) 

for all z G C \ Sfj,. 

Lemma 1. The vector-valued function is anti- analytic in C\S fl and has simple poles 
at points of S^. Moreover, it does not depend on the self-adjoint extension of A used to 
define ip(z). 

Proof. The first statement holds by rather obvious reasons, notice only that the poles 
of (ip(z),fi) coincide with those of ip(z). Let us pay attention to the second statement. 
Consider two different self-adjoint extensions Aft and AL. We have 

i/>(z) = (A s - z I) (Aft - ziy 1 % and i/>'(z) = (A\ - z l) (A\ - zl) i/> . 

Since both ip(z) and ip'(z) belong to Ker (A* — zl) and the dimension of this subspace is 
always equal to one, it follows that ip'(z) = g(z)tp(z) for every z G" Sp(A^) U Sp(Ai), where 
g(z) is a scalar function. Inserting this identity into (13. 2ft yields £(z) = £'( z )- D 

For every z G C \ we have the decomposition 7i = Ran (A — zl) + Span {/i} , in 
which case every element <p G Ti can be written as 

(p=[(p- tp{z)n\ + (p(z)n, 

where <p — <p(z)[i G Ran (A — zl). A simple computation shows that the non-orthogonal 
projection <p(z) is given by 

whenever z G C\S fl ; it is otherwise not defined. Indeed, because of the anti-linearity of the 
inner product in its first argument, the function <p(z) is analytic in <C\S fX and meromorphic 
in C for every ip G 7Y. We note that fl(z) = 1. 

Let us denote the linear map <p i— > (p(z) by $ M and the linear space of functions given by 
( 13.31) by $ M 7i. Since the operator A is simple, it follows that <1> M is injective and, therefore, 
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is an isomorphism from TC onto Qf/H |TTl Theorem 1.2.2]. Moreover, $ M transforms A into 
the multiplication operator on $^7^, that is, 

(Ay)(z) = (QpAQ- 1 !?) (z) = zfi(z), if G Dom(A). 

Proposition 1. Assume S^flR = 0. Let {x n } be the spectrum of any self- adjoint extension 
Aq of A. Then, for any analytic function f(z) that belongs to $^7^, we have 



X n 



T/ie convergence in flff.^p zs uniform over compact subsets of C\ 



(3.4) 



Proof. Fix some arbitrary self-adjoint extension Ajj of A. Take another self-adjoint ex- 
tension Ai ^ y4jj to define V ; ( 2 ); that is, ip(z) = (AJ — z /) (AJ — z/) ^q. Arrange the 
elements of Sp (A$) in a sequence {xnjneAf, where M is a countable indexing set, and 
let T] n be an eigenstate of A$ corresponding to x n , i.e., Afq n = x n r] n . Since A* D A$, 
it follows that ?] m G Ker (A* —x m I), where furthermore dim [Ker (A* — x m I)} = 1. On 
the other hand, since x m is not a pole of ip(z), i[>(x m ) is well defined and belongs also to 
Ker (A* — x m I). Therefore, up to a factor, r] m = £(x m ). 

Pick a sequence {Mk}kem of subsets of M such that C M^+i and UfcMfc = M. 
Consider any analytic function f(z) G $^7^. Since $ M is injective, there exists a unique 
ip ETC such that </?(z) = f(z). Clearly, 



neAf fc 



neMfc 
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The second factor in the r. h. s. of this inequality does not depend on z and obviously tends 
to zero as k — > oo. Since £(z) is continuous on C \ S^, the uniform convergence of (13.41) 
has been proven. □ 

Remark 2. Krein asserts that, for any operator in Symj^' 1 ^?^), one can always choose \i 
so that S M n E = pH Theorem 8]. 

Clearly, the sequence {x n } could be replaced by any other sequence {z n } of complex 
numbers for which {£,(z n )} is an orthogonal basis of 7i (for a related discussion, see [9]). In 
our case, the question of whether such a sequence exists or not is answered to the affirmative 
by invoking the self-adjoint extensions of the operator A. Notice also that, since for any real 
number x there is a self-adjoint extension of A containing x in its spectrum, it follows that 
every real point can be taken as a sampling point. In contrast with this, the construction 
of [8] has a fixed set of sampling points. 

Below we show that the interpolation formula (13. 41) is indeed a Lagrange interpolation 
series. 
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Proposition 2. Under the hypotheses of Proposition [7J there exists a complex function 
G(z), analytic in C \ (hence M.) and having simple zeroes at Sp(A^), such that 

/(*)= E 7 ^77 n /W, 

for every f(z) G $^7^. 



Proof. Given {x n } n€ M = Sp(Ajj), set ip(z) = (A$ — ZqI) (A$ — zl) 1 ip for some z such 
that Imzo^O and ip G Ker (A* - z I). Define 

G(S) : 



This function has simples zeroes at the poles of ip(z), that is, at points of the set {x n } n£ M- 
Also, £(z) = G(z)ip(z). Moreover, we can write 

tfj(z) = J!^L and G{z) = (z - x n )F{z), 

Z X n 

where rj(z) := [z — x n ) (A$ — z I) (A$ — zl) 1 tpo is analytic at z — x n , r)(x n ) ^ 0, and 
F{x n ) = G'(x n ). Thus, a straightforward computation shows that 

(£(*), £(x n )> G(z) (ri{z),r)(x n )) 



U(x n )\\ 2 (z-x n )G'(x n ) \\r)(x 



so it remains to verify that the last factor above equals one. By the Cauchy integral 
formula, we have 

V( x n) = — — <f dw 

= 77— f ( A t~ z I)(A^-wiy l tpodw 

^ 7r * J\w-x n \=e 

= -{A^-z I)(^-.j ^ (wl - A^y 1 dv?J ^ 

= ~ {X n - Zq) P n 4>0, 

where P n denotes the orthoprojector onto the eigenspace associated with x n . Therefore, 
(v(z),vM) = -(x n - Zq) (r](z),P n ip ) 

= -{X n - Zq) (Z - X n ) (P n {A^ - Z I) (A$ - zl)" 1 1p , i/)q) 

= \x n - z \ 2 (P n ipo,ipo) ■ 
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Finally, it is clear that (rj(x n ), r](x n )) = \x n — z \ 2 (P n ip ,ij; Q ) . □ 

Remark 3. Notice that the function G(z) is defined up to a constant factor. In particular, 
one may adjust it so that G'(xk) = 1, where Xk is a fixed eigenvalue of A$. Thus, a 
computation like the one in the proof above shows that 

G(z) = (z- Xk ) m ^ Xk » . (3.5) 
ll^)|| 2 

This identity may be useful in some applications; see Example 2 below. 



4 Spaces of analytic functions 

In this section we characterize the set of functions given by the mapping Note that in 
the general case Q^H is a space of analytic functions with simples poles in a subset of S^. 

Let 1Z C Ti be the linear space of elements (p for which (p(z) is analytic on 1R. As a 
consequence of [Til Corollary 1.2.1], it follows that 

/oo 
(f(x)rf(x)dm(x) (4.1) 
-oo 

for any ip,r] G 1Z and m(x) = (E x fi,fi), where E x is any spectral function of the operator 
A. That is, Q^IZ is a linear space of analytic functions in C \ 5* M such that their restriction 
to R belong to L 2 (IR, dm); in short, 

$^| K C L 2 (R,dm). 

Moreover, in this restricted sense $ M is an isometry from TZ into L 2 (IR, dm). 
The following theorem is due to Krein [191 Theorem 3]. 

Theorem 1 (Krein). For A G Syrn^ l \l~L), assume that 1Z = 7i. Consider a distribution 
function m(x) = {E x fi, fi), where E x is a spectral function of A. Then the map $ M generates 
a bijective isometry from 7i onto L 2 (R,dm) if and only if E x is orthogonal. 

This theorem deserves some comments. When E x is orthogonal it occurs that 

$ M ft| R = L 2 (R,dm) 

in the usual sense of equivalence classes. Thus, every function in L 2 (R,dm) is, up to 
a set of measure zero with respect to m(x), the restriction to 1 of a unique function 
that is the image under ^ of one and only one element belonging to 7i. Notice that 
any orthogonal spectral function of A corresponds to the spectral function of one of its 
self-adjoint extensions within 7i. Since these self-adjoint extensions have only discrete 
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spectrum, the inner product in L 2 (R, dm), with m(x) = {E x n, //), reduces to an expression 
like 

/■oo 

f(x)g{x)dm(x) = ^2c k f(x k )g(x k ) 

k 

whenever E x is orthogonal. That is, the equivalence classes in these spaces are quite broad. 

The following corollary is partly a straightforward consequence of Theorem [TJ Notice 
that 5 M nR = implies TZ = H. 

Corollary 1. Let A G Sym^ (TC) and choose a gauge /i for this operator. Assume that 
n R = 0. Let E x be one of its orthogonal spectral functions. Then the linear space of 
functions TC^ := Q^rt, equipped with the inner product 



(/(•)> #(')) : = / f(x)g(x)dm(x) where m{x) = (E x fi, //) , (4.2) 

«/ — oo 

zs a reproducing kernel Hilbert space, with reproducing kernel k(z,w) := £(u>)). 
Proof. We only verify the last statement. Given /(z) = (£(z),<£>) G 7-^, we have 



(K, w )J(-))= / k(x,w)f(x)dm(x) = (£(w),<p) = f(w), 

J — oo 

where the second equality follows from ( 14. ip . □ 

Notice that, once the operator A is given, the linear space 7i M depends only on the 
choice of gauge \i. By Corollary [I], for those gauges that obeys fl R = 0, 7Y M may be 
endowed with different Hilbert space structures, one for each orthogonal spectral function 
of the operator A. By ( 14. ip . all these Hilbert spaces are however isometrically equivalent. 

Irrespective of anyone of these Hilbert space structures, TC^ possesses the following 
properties. 

Proposition 3. Suppose \i such that fl R = and consider 7i M equipped with an inner 
product of the form fl^.jgp . 

(i) Let w be a non real zero of f \z) G H^. Then the function g(z) := f(z)(z—w)(z — w)^ 1 
is also in Ti^ and \\g{-)\\ = ||/(-)ll- 

(ii) The evaluation functional, defined by /(•) t— > f(z), is continuous. 

Proof, (i) Suppose that w is a non real zero of f(z). Since f(z) = (£(z), ip) for some (p <EH, 
it follows that ip is orthogonal to ip{w) and therefore ip G Ran (A — wl). Note that 

f(z) = (fa), (A - wl) (A - wl)' 1 cp) = {z- w) (C{z), (A - wl)' 1 ip) . 
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In these computations we have used that £(z) G Ker (A* — z). Moreover, 

(£(z), (A - wl) (A - wiy 1 <p) = (z- w) (£(*), (A - wiy 1 ip) = —/(*). 

Then /(z) (2 — w) (z — G 7"^. The equality of norms follows from the fact that the 

Cayley transform — wl) (A$ — wl)^ 1 is an isometry. 

(ii) Let f(z),g(z) G Hp. Then f(z) = (£(z),ip) and g(z) = (£(z),r)) for some <p,r}€H, 
and furthermore 

\f(w)- 9 (w)\ = \(t{w),<p-r,)\ < umw y-v\\ = wtww ii/(o -0(011 • 

In other words, this result follows from the fact that H^ is a reproducing kernel Hilbert 
space. □ 

Definition 2. An operator A G Sym^' 1 ' ) (7i) is called entire if there exists a so-called entire 
gauge ji <EH such that <p(z) is an entire function for every ip G H. Equivalently, A is entire 
if H = Ran (A — zl) + Span {/x} for all z G C. 

Notice that if A is entire, £(z) is a vector- valued entire function and is a Hilbert 
space of entire functions. The spaces of reconstructible functions treated in [8] corresponds 
to a special situation within this particular case of our approach. 

Definition 3. A Hilbert space of entire functions is called a de Branges space if, for every 
f(z) in that space, the following conditions holds: 

(i) for every non real zero w of f(z), the function f(z)(z — w)(z — w)' 1 belongs to the 
Hilbert space and has the same norm as f(z), 

(ii) the function f*(z) := f(z) belongs to the Hilbert space and also has the same norm 
as f(z); 

and furthermore, 

(iii) for every w:Imw/0, the linear functional f(-) 1— > f(w) is continuous. 

There is an extensive literature concerning the properties of de Branges spaces. We 
refer to [3] for more details. 

It can be shown that H^ is a de Branges space for certain choices of the entire gauge 
Ii. There are some evidence indicating that Krein noticed this fact [HI pag. 209]. Also, 
some hints supporting this assertion has been given by de Branges himself [5]. We however 
could not find any formal proof of this statement. Thus, for the sake of completeness and 
for the lack of a proper reference, we provide a proof below. 
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Remark 4. As a consequence of [TIJ Lemma 2.7.1], given any self-adjoint extension A$ of 
an operator A G Sym^' 1 ' ) (H), one can always find a complex conjugation C for which A$ is 
real. If follows from the proof of the cited lemma that Cip(z) = ip(z) when ip(z) is written 
in terms of the real self-adjoint extension. Moreover, by [TTl Theorem 2.7.1], the operator 
A is also real with respect to C . Since by [23l Corollary 2.5] all the self-adjoint extensions 
of A are real, it follows that Cip(z) = ip(z) for every realization of ip(z). 

If furthermore A is entire, then an entire gauge fi may be chosen real, that is, C fi = fi 
(see [m Theorem 1] and also [TY], Sec. 2.7.7]). 

Proposition 4. Assume that an entire operator A is real with respect to some complex 
conjugation C and let fi be a real entire gauge. Then the associated Hilbert space TC^ is a 
de Branges space. 

Proof. In view of Proposition [21 we only have to verify (ii). By Remark H] we know that 
Ci/j(z) = ip(z) thence C£(z) = £(z). Now consider any f(z) = (£(z),(p). Clearly f*(z) : = 
(£(z),Ccp) also belongs to Ti.„. Furthermore, 



f*(z) = (t{z),C<p) = (C£(z), V ) = (t®, V ) = f(z). 

Since C is an isometry, the equality of norms follows. □ 

Notice that we only have used the simultaneous reality of the entire operator and its 
entire gauge \i in showing that TC^ obeys (ii) of Definition [3j Indeed, this condition is also 
necessary. 

Proposition 5. IfH^ is a de Branges space there is a complex conjugation C with respect 
to which both A and fi are real. 



Proof. Let C : — > be defined by (Cf)(z) = f(z), for every f(z) G H^. Clearly, C is 
a complex conjugation. Moreover, the multiplication operator A, defined on the maximal 
domain in TC^ by (Af)(z) = zf(z), is real with respect to C. Now define C := 
By construction, C is a complex conjugation in H^. Since (CJt)(z) = 1 = Ji(z), it follows 
that fi is real with respect to C. Finally, it is not difficult to see that A = and 
therefore A is real with respect to C . □ 

An alternative definition of a de Branges space is given in terms of an arbitrary entire 
function e(z). The de Branges space TC(e) associated to e(z) is 



H{e) = if{z) G Hol(C) : f(z)/e(z) G A/", f{z)/e{z) G A", J' 



\f(x)/e(x)\ 2 dx < oo 



where M denotes the class of analytic functions of bounded type and non-positive mean 
type in the upper half plane [H Chap. 1]. From this definition, it follows that H{e) is a 
reproducing kernel Hilbert space, whose reproducing kernel is expressed in terms of the 
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function e(z). As shown in [U Chap. 2] (see also [T2J Sec. 5]) a Hilbert space of entire 
functions 7i that obeys Definition [3] is unitarily equivalent to a de Branges space Tl(e) 
with 

e{z) = ij — \— — r(u% - z)k(w , z), 

V k(w ,w ) lm(w ) 

where k(w, z) is the reproducing kernel of 7i and Wq is any non-real complex number such 
that k(wo, wo) > 0. 

As discussed in [U [12], the multiplication operator A with maximal domain in TC(e) is 
a closed, symmetric operator with deficiency indices (1,1) and domain of codimension 
or 1. Since in our case the multiplication operator is unitarily equivalent to an operator 
belonging to Sym^'^(H), the codimension of Dom(A) is necessarily equal to 0, that is, A 
is densely defined. The set of self-adjoint extensions of A are in one-one correspondence 
with the set of entire functions s t (z) := — sinta(z) + costb(z), t G [0, 7r), where a(z) and 
b(z) are defined by the identity e(z) = a(z) + ib(z). In terms of s t (z), we have 

Dom(A) = \g(z) = f( Z ) ~ s t (z) f (w ) g ^ ; _^ 1 

A^(z) = zg(z) + f(w )s t (z). 

For details, see [121 Sec. 6]. Notice that the definition of Dom(,4 t ) does not depend on the 
choice of wq, as one can verify by resorting to the first resolvent identity. 



5 Examples 

In this section we give two elementary illustrations of the method developed above. 
1. Consider the semi- infinite Jacobi matrix 



( ql 


h 








bi 




b 2 








b 2 




h 








h 













(5.1) 



•7 



where b^ > and e M for fceE Fix an orthonormal basis {5fc}fc S N in 7i. Let J be the 
operator in Ti whose matrix representation with respect to {^j^gN is ( 15.11) . Thus, J is the 
minimal closed operator satisfying 



(5 n , J5 n ) = q n , (S n+ i, J5 n ) = (5 n , J5 n+1 ) 



VneN 



(Consult [2 Sec. 47] for a discussion on matrix representation of unbounded symmetric 
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operators.) It is well known that J may have either deficiency indices (1,1) or (0,0) [TJ 
Chap. 4, Sec. 1.2]. A classical result is that if J has deficiency indices (1,1), then the 
orthogonal polynomials of the first kind Pk{z) associated with (15. ip are such that 



Y,\Pk(z)\ 2 < oo 



fc=0 



uniformly in any compact domain of the complex plane [TJ Theorem 1.3.2] . Therefore, for 
any z G C, n(z) = Y^T=i Pk-i{z)$k is i n H*. By construction, ir(z) is in the one- dimensional 
space Ker(J* — zT). It is also known that, when the deficiency indices are (1,1), J is an 
entire operator and Si is an entire gauge for J pj], Sec. 3.1.1 and Theorem 3.1.2]. 

Let us find £(z) for the operator J. Taking into account (13.21) . (5i,£(z)) = 1 and 
(5i,n(z)) = 1 for all z G C. Then, since both 7r(z) and £(z) are in Ker(J* — zl) and Si 
is entire, ir(z) = £(z) for all z G C. Thus, for any ip in Ti, <p = J2'k'=i i PkSk, the function 
tp(z) G H Sl is given by 

oo 

(p{z) := (n(z),<p) = ^P k -i{z)(p k , zEC. 

k=l 



Clearly, if (p(z) G 7i,s 1 , <fi(z) G TCs 1 - Whence, in virtue of Proposition [3] and, our 
space TCs 1 is a de Branges space and then, by Proposition [5l Si is real with respect to 
C = (C is the conjugation in Ti,s 1 given in the proof of Proposition [5]). 

Formula (13.41) is written in this 



x„gSp( J» 



^ "^n)|| 2 ^ 



XnGSpl Jt 



where Jj is a certain self-adjoint extension of J. 

In a different setting, sampling formulae obtained on the basis of Jacobi operators have 
been studied before [6j [7j. We remark that in [6] the interpolation formulae differ from 
the ones obtained by us. In [7], Jacobi operators are treated without using M.G. Krein's 
theory of entire operators. 

2. The entire operator used here has been taken from [TT] and is a particular case of an 
example given by Krein in [2Tj . 

Consider a non- decreasing bounded function s(t) such that 

s(-oo) = and s(t - 0) = s(t) . 
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Fix a function defined for any x in the real interval (—a, a) by 



F(x) := / e ixt ds{t) . 

J — oo 

In the linear space L of continuous functions in [0, a) vanishing in some left neighborhood 
of a, we define a sesquilinear form as follows 

pa pa 

(g,f):= / F(x-t)f(x)yjtjdxdt. (5.2) 
Jo Jo 

This form is a quasi-scalar product, i. e., the existence of elements / in C such that / ^ 
and nevertheless (/, /) = is not excluded. 

Denote by V the set of continuously differentiable functions / £ C such that /(0) = 
and define in T> the differential operator A by the rule Af := if. It is not difficult to show 
that (g,Af) = (Ag,f) and V is quasi-dense in C Now, proceeding as in [11], Sec. 2.8.2], 
one defines the space C as follows 

C = C\0, = {/e£:(/,/) = 0}. 
In we define an inner product by 

(V,<f):=(g,f), (5.3) 

where <^ and 77 are equivalence classes containing / and g, respectively. Let TC be the 
completion of £ and consider in it the operator A such that, for the equivalence class ip 
containing / £ T>, Atp is the equivalence class containing Af. It can be shown that A is 
an entire operator and that 

(p(z) = (£(z),<p)= [ a e izt f(t)dt, 
Jo 

where / £ (p. This identity, together with (15.21) and ( 15.31) . determines £(2) completely [HI 
Sec. 3.2.2]. 

Notice that, in this example, the entire gauge associated with £(z) remains unknown. 
This is not however an issue since the sampling kernel can be computed anyway by resorting 
to expression (13. 5p . 

Acknowledgments. We express our gratitude to M. Ballesteros for drawing our attention 
to the results of [13] . We also thank A. Minkin for useful remarks that led to an improved 
version of this work. 
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